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Abstract. Standard perturbative calculations lead to pathologically large NLO corrections to 
low- a; Bj evolution equations like BFKL and BK. Using a more refined treatment of kinematics 
in mixed-space, relevant when gluon saturation sets on, one obtains an improved version of the 
BK equation, resumming to all orders the most severe of those large higher order corrections. 



1. Introduction 

Hadronic collisions at very high energy are involving partons with very small momentum fraction 
in the hadronic wave- functions. Due to the high occupancy of those wee partons (mostly 
gluons), the phenomenon of gluon saturation occurs: multiple parton scattering is typical, and 
accompanied by strong color coherence effects. Hence, the collinear factorization (and other 
standard perturbative QCD formalisms) which involves only one parton from each colliding 
hadron does not capture the typical physics of collisions in the high-energy limit, with only a 
semi-hard momentum transfer. 

Instead, the coherent multiple scattering effects are taken into account straightforwardly 
when describing the wee gluons inside each highly boosted hadron as a semi-classical gluon field 
(see Ref. pQ and references therein). The main formalism based on this idea is the Color-Glass- 
Condensate effective theory (CGC). Each ultra-relativistic nucleus is described by a random 
classical shockwave gluon field with a classical statistical distribution, and QCD quantum 
corrections are resummed within leading logarithmic accuracy (LL) at small momentum fraction 
by the JIMWLK evolution of the gluon field distribution. The JIMWLK functional equation 
can also be written as Balitsky's infinite hierarchy of equations. 

Inclusive enough observables, like DIS structure functions at low XBj or single inclusive 
particle production at forward rapidity in pA collisions, can be expressed in terms of the 
scattering amplitude of a color dipole on the gluon field of the target. For that object, the 
Balitsky-Kovchegov (BK) equation [2j El E] gives a safe approximation of the full JIMWLK 
equation. Adding running coupling effects to the BK equation [5j [6] leads to a successful 
phenomenological description of DIS data at HERA 0[8] within the CGC, as well as of forward 
particle production at RHIC [9]. 

In the recent years, the calculation of NLO corrections in that framework has been a hot topic. 
After the running coupling corrections [5j[6], the full set of NLO corrections to the BK equation 
have been calculated |10j . Later, the calculation of the NLO corrections has been performed for 



the impact factor or coefficient function both in the case of DIS structure functions |11] [T2] and 
of forward hadron production in pA collisions [13] . 

Unfortunately, those NLO calculations cannot be used, in the form in which they are now 
available, to perform phenomenological studies at full NLO accuracy. Indeed, the BK equation 
at NLO suffers from the same problem as its linear version, the BFKL equation at NLO [14] 115] : 
some of the NLO corrections are pathologically large and lead an instability of the solutions. 
This signals a breakdown of the perturbative expansion as done usually in the Regge limit. 
Those large NLO corrections are due to the inability of the standard perturbative expansion 
in the Regge limit to provide results matching smoothly with DGLAP physics in the collinear 
and in the anticollinear regimes [TB]. Hence, the large higher order corrections to the BFKL 
and BK equations can be resummed to all orders by performing an appropriate matching with 
the DGLAP equation at LO (or beyond) in the collinear and in the anticollinear regimes |16j . 
That program has been completed for the BFKL equation both in momentum space [T7] and 
in Mellin space |18| . However, the BK equation is more naturally written in mixed space. The 
generalization of that resummation to the case of the BK equation requires a significant effort 
mostly due to the translation to mixed space, and to a lesser extent due to the nonlinearity of 
the BK equation. 

Among the large higher order corrections to be resummed, the ones of purely kinematical 
origin are the most severe, but also the easiest to deal with [16J. In section[2]of this contribution, 
the NLO impact factors for DIS [12] are analysed to understand why such kinematical issues 
arise. And in the section [31 an improved version of the BK equation at LO is proposed, which 
realizes the resummation of those large kinematical higher order corrections. It corresponds to 
the mixed spac^l analog of the kinematical constraint |19l 1201 I2T] in momentum space. It also 
represents a first step towards a full resummation providing a fully stable and reliable version 
of the BK equation at NLO. 



2. Diagnosing kinematical issues from the explicit NLO impact factors for DIS 

The DIS structure functions are linear combinations of the total cross sections for the scattering 
of a transverse or longitudinal virtual photon off the target, which at strict NLO accuracy in 
the CGC can be written as [T2] (see also [TT] ) 
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where Sij is the S-matrix for the scattering of a fundamental color dipole with transverse 
positions x, and Xj off a gluon shockwave, (. . .}o is the statistical average over the target's 
gluon field with no LL quantum corrections included, and Xij = |xj — Xj\. The LO impact 
factors Z^i have been known for a long time [22] I23j . whereas the NLO ones have been 

calculated in Refs. |11[ 112]. The integral over the photon's momentum fraction z 2 = k^/q^ 
carried by the radiated gluon is logarithmically divergent for z 2 — > 0, but in that limit 

2 

l^°(x 0l x 1 ,x 2 ,z 1 ,z 2 = 0) = J^ r TL° L (x x,z 1 ). (2) 

x 02 x 2l 



1 In mixed space, the kinematics of partons is described by their light-cone momentum k + and their transverse 
position x. 



Together with an appropriate factorization scheme (including for example a cut-off in k + ), the 
BK equation 
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allows to resum those small 2:2 LL contributions. In that case, one should use in the first line 
of the expression (P) the dipole S-matrix (Soi) Y + evolved with the BK equation ([3|) over a 

range = log(fcj /fc+ in ). fc^ in is the typical k + scale set by the target and kf an appropriate 
factorization scale in k + , such as = Z\{1 — z\)q + . 

The LO (resp. NLO) impact factor (resp. T^jf) contains a factor which suppresses 
exponentially the large values of Q 2 Xf (resp. Q 2 X 2 ), where 

X 2 = z\ (l-zi)xo! and X 2 = z 1 (1- z 1 - z 2 ) x^ + z 2 (I- z ± - z 2 ) xj 2 + z 2 z 1 x^ . (4) 

As argued in Ref. [12] the variables Q 2 X 2 and Q 2 X 2 are the ratios of the formation time 
of the quark-antiquark or quark-antiquark-gluon Fock components of the photon, resolved by 
interaction with the target, over the lifetime of the virtual photon. Hence, the interpretation of 
that exponential suppression is very clear: a Fock state which has not enough time to be formed 
as fluctuation of the virtual photon within the lifetime of the latter cannot give a non-negligible 
contribution to the DIS cross sections. 

The standard treatment of low z 2 LL with the BK equation discussed previously requires 
to approximate X 2 by X 2 , in order to obtain the factorization ([2]) of Z^Y£°. Although exact 
at z 2 = 0, the approximation X 2 ~ X 2 is not generically correct at small but finite z 2 : it is 
wrong when the gluon is emitted at a so distant transverse position x 2 that zi(l — 2i)zqi "C 
22^02 — z 2 x i 2 - In that regime, not only the nice feature of suppression of Fock states too long 
to form is spoiled by the standard subtraction of LL, but also the term used to subtract the LL 
contributions from the NLO term in the expression ([1]) is parametrically larger than both the 
unsubtracted NLO term and the LO term, which signals a breakdown of this formalism. 

Evolution equations like BK and BFKL can be derived from the knowledge of the photon 
impact factor at arbitrary order but restricted to the case of softer and softer gluons emitted 
successively \2A\ [3] . Usually, all the transverse scales are assumed to be of the same order in that 
context. This assumption (also used in other derivations of those equations) is not completely 
self-consistent due to the unrestricted integration over transverse momentum or position in the 
evolution kernel. The resulting issue is essentially the same as found in the study of the NLO 
photon impact factor: in the parton cascades resummed by the LO BFKL and BK equations, 
the softer and softer gluons are not always correctly ordered in formation time. The most 
pathological higher order corrections to the BFKL and BK kernel are then induced by that little 
inconsistency at LO. 

3. Improving the treatment of kinematics in the BK equation 

According to our previous discussion, the standard BK equation at LO ([3]) includes, at large 
x 2 , unphysical contributions from gluons which should not have time to be formed. The first 
step to cure this problem is to modify the probability density of soft real gluon emission by a 
color dipole by forbidding emissions with Zi(l— z^x^ <C z 2 Xq 2 — z 2%i 2 - This is the mixed-space 
analog of the kinematical (a.k.a. consistency) constraint of Refs. [191 1201 f2Tj . In mixed-space, 
this restriction has been first proposed in Ref. [25], where it was inferred from the structure of 
n-gluons MHV amplitudes. The resummation scheme proposed here differs however from the 
one in Ref. [25] in several aspects, most notably in the treatment of virtual corrections. 

The choice done in Ref. [26J is to organize the virtual corrections in such a way that the 
probabilistic interpretation of the dipole cascade |24] is maintained at each order once the 



resummation is done. Then, the modified virtual corrections can be calculated unambiguously 
from the modified real gluon emission probability. Performing this task and writing the result 
as an evolution equation, one gets the improved BK equation with kinematical constraint |26j 
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A convenient definition for the shift A012 (but not unique, due to some resummation scheme 
ambiguity) is 

A 012 = max|o, log ^ ^fe^l H , (6) 

Due to the theta function in ([5]), the phase space for gluon emission at large X2 is severely 
restricted at small Y + , but progressively opens up in the course of the Y + evolution. The shift 
of Y + in the real term only also contributes to slow down the evolution with respect to the 
standard BK equation Q. 

The largest NLO corrections [10] to the BK equation are indeed resummed into the improved 
LO equation ([5]). However, a full resummation of all pathological NLO corrections as in Refs. 
[TTJ, [18] requires further work. The equation ([5]) also allows to subtract the LL contributions 
from the NLO impact factors ([T]) in a correct way in all the phase-space, by contrast to the 
equation 

Implementing the kinematical improvement (0) together with running coupling corrections 
[5] should lead to solutions more naturally in good phenomenological agreement with the DIS 
data, due to a slower Y + evolution [HE]. 
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